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Abstract. We study a simple stochastic differential equation that models the dispersion of close heavy 
particles moving in a turbulent flow. In one and two dimensions, the model is closely related to the one- 
dimensional stationary Schrodinger equation in a random <5-correlatcd potential. The ergodic properties of 
£N| i the dispersion process are investigated by proving that its generator is hypoelliptic and using control theory. 
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Cn ' 1. INTRODUCTION 

Transport by turbulent flows belongs to phenomena whose understanding is both important for practical 
*£h applications and abounds in intellectual challenges. Unlike the reputedly difficult problem of turbulence 

per se, turbulent transport allows simple modeling that accounts, at least qualitatively, for many of its 
■ observable features. The simplest of such models study transport properties of synthetic random velocity 
fields with presupposed distributions that only vaguely render the statistics of realistic turbulent velocities. 
The advection by velocity fields of quantities like temperature or tracer density may be derived from the 
dynamics of the Lagrangian trajectories of fluid elements. In synthetic velocity ensembles, such dynamics 
is described by a random dynamical system. One of the best studied schemes of this type is the so called 
Kraichnan model based on a Gaussian ensemble of velocities decorrelated in time but with long-range 
spatial correlations [HUH] • In this case, the random dynamical system that describes the Lagrangian flow is 
qq given by stochastic differential equations (SDE's). It was successfully studied with the standard tools of the 

£ — . • theory of random dynamical systems, but it also led to non-trivial extensions of that theory [T2irT9"ll20] . 

The problem of turbulent transport of matter composed of small but heavy particles (like water droplets in 
ON ' turbulent atmosphere) may be also studied by modeling turbulent velocities by a random synthetic ensemble, 
but it requires a modification of the previous approach. The reason is that heavy particles do not follow 
Lagrangian trajectories due to their inertia. On the other hand, the assumptions of the time decorrelation 
of random velocities may be more realistic for inertial particles on scales where the typical relaxation time 
of particle trajectories (called the Stokes time) is much longer than the typical correlation time of fluid 
velocities. There have been a number of papers that pursued the study of dynamics of inertial particles 
with various simplifying assumptions, see e.g. [Ttt71I^1 [T71 122112511291133] ■ The primary focus of those studies, 
combining analytical and numerical approaches, was the phenomenon of intermittent clustering of inertial 
particles transported by turbulent flow. A good understanding of that phenomenon is of crucial importance 
for practical applications. 

The aim of the present article is to show that the simplest among the models of inertial particles dynamics 
are amenable to rigorous mathematical analysis. More concretely, we study the SDE's that describe the pair 
dispersion of close inertial particles in shortly correlated moderately turbulent homogeneous and isotropic 
G?-dimensional velocity fields (not necessarily compressible). Such models were discussed in some detail 
in [2j[T7l|23 ( 24, 29, 33 . In particular, it was noted in [33] that the d = 1 version of the model is closely 
related to the one-dimensional stationary Schrodinger equation with (^-correlated potential studied already 
in the sixties of the last century [2] as a model for Anderson localization. As was stressed in [T7] . the 
d = 2 model for the inertial particle dispersion is also related to the one-dimensional stationary Schrodinger 
equation, but this time with <5-correlated complex potential. The models for dispersion were used to extract 
information about the (top) Lyapunov exponent for the inertial particles which is a rough measure of the 
tendency of particles to separate or to cluster [21133] . The numerical calculations of the Lyapunov exponents 
in two or more dimensional models of particle dispersion presumed certain ergodic properties that seemed 
consistent with results of the simulations. Here, we shall establish those properties rigorously by showing 
the hypoellipticity of the generator of the Markov process solving the corresponding SDE and by proving 
the irreducibility of the process with the help of control theory. For a quick introduction to such, by now 
standard, methods, we refer the reader to [TlJ[30]. More information about the ergodic theory of Markov 
processes may be found in the treatise [25]. 

l 



2 



KRZYSZTOF GAWEDZKI, DAVID P. HERZOG, AND JAN WEHR 



The paper is organized as follows. In Sec. [21 we present the SDE modeling the inertial particle dispersion. 
In Sec.EJ we recall its relation to models of one-dimensional Anderson localization. Sec.|4] establishes the 
hypoelliptic properties of the generator of the dispersion process. In Sec.El we introduce the (real- projective 
version of the dispersion process whose compact space of states may be identified with the (2d— 1 )-dimensional 
sphere S 2 ^ 1 . Sec. [5] is devoted to proving that the dispersion process is controllable. Together with the 
hypoelliptic properties of the generator, this implies that the projectivized version of the process has a 
unique invariant probability measure with a smooth strictly positive density. The analytic expression for such 
a measure may be written down explicitly in d — 1 but not in higher dimensions. The smoothness and strict 
positivity of its density provides, however, in conjunction with the isotropy assumption, valuable information 
about the equal-time statistics of the projectivized dispersion. The isotropy permits to project further the 
projectivized dispersion to the quotient space 5 2d_1 /50(rf). For d = 2, this space may be identified with 
with the complex projective space PC 1 = CUoo and the projected process with the complex-pro jectivized 
dispersion. In fact, the point at infinity may be dropped from PC 1 since the complex-projectivized dispersion 
process stays in C with probability one. For d > 3, the quotient space S 2d ~ 1 /SO(d) is not smooth but has 
an open dense subset that may be identified with the complex upper-half-plane that the projected process 
never leaves. These non-explosive behaviors are established in Sec. [8] by constructing a Lyapunov function 
with appropriate properties. Finally, in Sec. [5] and Appendix [Bl we demonstrate how the established ergodic 
properties of the projectivized dispersion process lead to the formulae for the top Lyapunov exponent for 
inertial particles that were used in the physical literature. Appendix [Al derives a formula, used in the main 
text, that expresses the S'0(2d)-invariant measure on S 26 " -1 in terms of SO(d) invariants. 

Acknowledgements. J.W. acknowledges the support from C.N.R.S. for an extended visit at ENS-Lyon 
during which this work was started. D.P.H. and J.W. acknowledge partial support under NSF Grant DMS 
0623941. 



2. BASIC EQUATIONS 

The motion in a turbulent flow of a small body of large density, called below an inertial particle, is well 
described by the equation 



r = -J(r-u(i,r)), (2.1) 

where r(t) is the position of the particle at time t and u(t,r) is the fluid velocity field. Relation (|2.ip is 
the Newton equation with the particle acceleration determined by a viscous friction force proportional to 
the relative velocity of the particle with respect to the fluid. Constant t is the Stokes time. Much of the 
characteristic features of the distribution of non-interacting inertial particles moving in the flow according 
to Eqs. (|2.ip is determined by the dynamics of the separation 8r(t) = p(t), called particle dispersion, 
of very close trajectories. In a moderately turbulent flow, the particle dispersion evolves according to the 
linearized equation: 

p = -±(p-(p-V)u(t,r(t))) (2.2) 

or, in the first-order form: 

P = 7X, X = ~X + (P- V)u(t,r(t)). (2.3) 

For sufficiently heavy particles, the correlation time of (Vtt)(i, r(t)) is short with respect to the Stokes time 
r and one may set in good approximation [2] 

V jU l {t,r{t))dt = dS){t) (2.4) 

where dS(t) is a matrix- valued white noise with the isotropic covariance 

(dS l j(t) dSf(t')) = D% S(t - t') dtdt' , Djf =AS ik 8 J i + B{5}5\ + <5j<$) . (2.5) 

Positivity of the covariance requires that 

A>\B\, A + (d + l)B>0. (2.6) 

Incompressibility implies that A+ (d+l)B — 0, but we shall not impose it, in general. We shall only assume 
that A + 2B > for d = 1 and that A > for d > 2. After the substitution of (f2H) . Eq. (j2~2"j) becomes 
the linear SDE 

1 . . 1 dS(t) , , 

P = --P+ ~—P (2-7) 
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that may be written in the first order form in a more standard notation employing differentials as 

' P \ _ ( o \dt \ ( 

We shall interpret the latter SDE using the Ito convention, but the Stratonovich convention would lead to 
the same process. The solution of Eq. (|2.8[) exists with probability 1 for all times and has the form 



(2.8) 



X (t)J ~ XeXP Ly \dS{s) -\ds)\ { X (0))' 



(2.9) 

where the time ordered exponential may be defined as the sum of its Wiener chaos decomposition 



T exp 



„-.„,(; a) 



±ds 
dS{s) -±ds 



E 



n ~° 0<a 1 <—<s n <t 





ds 1 ...ds n (2.10) 



0\ 4 ) ( 0\ 31 (o 4) 
v dS(s„) oj 6 "' \dS( 3l ) Oj 6 

that converges in the L 2 -norm for functionals of the white noise dS(t). The resulting stochastic process 
(p(t), x(t)) = P(t) i s Markov and has generator 

1 = 7(rv p - r vj + ^ P V^|v x ,v xfc . (2.11) 

i,j,k,l 

In other words, for smooth functions /, 

|(/(P(*))) = (2- 12 ) 

for ( • • • ) denoting the expectation. 

For the process p(i) given by Eq. ([2l9]l . p(i) = for all < > if p(0) = 0. On the other hand, if p(0) 7^ 
then p(t) 7^ with probability 1 for all t > so that we may restrict the space of states of the Markov 
process p(t) to R 2d \ {0} = R^ . 

3. RELATION TO ONE-DIMENSIONAL LOCALIZATION 

In the lowest dimensions, there is a simple relation between the stochastic process p(t) = (p(t),x{t)) 
constructed above and simple models of Anderson localization in one space dimension. Let us set 

4>(t) = e^p(t) (3.1) 

exponentially blowing up the long-time values of pit). Eq. ()2.7[) implies that 

"^ + 7^ = -^, (3.2) 

or, in the first order form, 

S) = {ds(t)°+±dt T o) (t) ' (3 ' 3) 

Similarly as before, the above SDE defines a Markov process. Clearly, 

mt),m = ^{p{t), X {t)+\p{t)). (3.4) 

Viewing t as the one-dimensional spatial coordinate, Eq. (|3.2[) takes the form of the vector-like stationary 
Schrodinger equation 

- J^/> + V(t)ip = Eip, (3.5) 

where V(t) — - dS ^' plays the role of the random matrix- valued white- noise potential and E = —jipi of the 
(negative) energy. In particular, in d = 1, ip(t) is a real scalar function and so is the (5-correlated potential 

V(t) = I^AT^^P, (3.6) 
r at 
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where f3(t) is the Brownian motion. The scalar version of Eq. Q3.5P was studied in [14] as a model of one- 
dimensional Anderson localization, see also [21]. In d = 2, interpreting ij) as a complex number tp 1 + itp 2 , 
one may replace the matrix valued (^-correlated potential V(t) in the SDE (|3.5|) with the complex valued 
one 

where /3 1 (i), /? 2 (i) are two independent Brownian motions (the two realizations of V{t) lead to the Markov 
processes with the same generator and, consequently, with the same law). Consequently, as stressed in [17) . 
Eq. p.5p in d = 2 may be viewed as a model of localization for a one-dimensional non-hermitian random 
Schodinger operator of the type not studied before. 

4. HYPOELLIPTIC PROPERTIES OF THE GENERATOR 

The generator (12. lip of the particle dispersion process p(t) = (p(t),x(t)) has certain non-degeneracy prop- 
erties which imply smoothness of the transition probabilities. Let us start with the following fact about the 
covariance (|2.5I) of the matrix- valued white noise dS(t) : 

Lemma 4.1. We have 

d 

Df( = (ES'S™ + F5 vm S jn + G5* n S™)(E5fS™ + F5 krn S ln + G5*5r) (4.1) 

m,n— 1 

for 



2 

Proof. The right hand side of Eq. (|4.1[) is 



E = d~ l ( - VA + B + y/A+{d+\)B) , 
F = ^{VA + B + VA-B), 

G = UVA + B- VA^B). (4.2) 



{E 2 d + 2EF + 2EG)8)8 k + (F 2 + G 2 )S ik S :j i + 2FG8\5 k . (4.3) 



Hence, in order to satisfy Eq. (|4.ip . we must have 

A = F 2 +G 2 , B = E 2 d + 2E{F + G)=2FG. (4.4) 
The assumed values of E, F, G solve these equations. □ 
Define the vector fields 

Ac = |(X-V P - X -V X ), 

Y™ = -Y,(^{E6)5^ + FS im 5 jn + G5i5Y)V pi 

i,3 

+ E(x j + ^{E^Sn + FS* m 8 jn + Gtf n fi?) V x< = r[A , A™] , 

Z? = - ^( V + fj){E8}8™ + F8™8 ]n + G8i5f) V p * 

id 

+ + ^M-C + + V x * = r[A , r„ m ] ■ (4-5) 



Using Eq. (|4.ip . one infers that Eq. (|2.11l) giving the generator L may be rewritten in the form: 

l = xo + !E(*™) 2 - ( 4 - 6 ) 

Remark 4.2. The process p(t) may be equivalently obtained from the SDE 

d 

dp = X (p)dt+ J2 Xn(p)d/3l(t), (4.7) 

m,n— 1 

where are independent Brownian motions. Here, we adopt the Stratonovich convention and hence the 

generator corresponding to that equation has the form (|4.6p so that the latter SDE leads to a process with 
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the same law as pit). The convention of the stochastic integral we choose, however, is insignificant as the 
process p[t) is of bounded variation. 

In order to establish hypoelliptic properties of L, we shall use the following non-degeneracy relation 
satisfied by the vector fields X™, Y™ and Z™ : 

Proposition 4.3. Suppose that p = (p, X ) ^ 0. Then the vectors 

X™(p), Y™(p), Z™(p) with m,n=l,...,d (4.8) 

span the 2d-dimensional space. 

Proof. First suppose that p = so that x 7^ 0. We have 

Yn (0,X) = E ( ES i 6 n + FSim ^n + GSffi) V x « . (4.9) 

Let e R d . Then 

E(«XY"+/5XT)C(0,X) = [E(ax 2 +f3x-<P) + Fax 2 + G(a X 2 +Px-^)}x- V x 

m,n 

+ F/3 X 2 0-V X . (4.10) 

Setting 



(E + F + G)F (x 2 ) 2 ' ^X 2 
(note that F > and E + F + G > 0), we obtain 

E(«xV n + /5xV n )^ m (0,x) = 4>-V x . (4.12) 



Hence the vector • V x is in the span of (|4.8|) for arbitrary cf>. We have still to show that an arbitrary 
vector <j ■ V p is in that span. To this aim note that 

Z™(0, X ) = -2j2x J (ES}6™ + F5* m S jn + GS^Sf) V p , + ^-^r . (4.13) 

i,j i 

Proceeding as before, we show that an appropriate combination of -Z™(0, x) gives the vector 

<x.V p + ]T(-rV x! (4.14) 

i 

from which the term XX-O^x* ma y De removed by subtracting an appropriate combination of Y™(Q, X ). 
That ends the proof of the claim of Proposition 14.31 for p = 0. 

Suppose now that p ^ 0. Proceeding as before, we see that arbitrary vector • V x may be obtained 
by taking an appropriate combination of the vectors X™(p, %). Similarly, arbitrary vector cr ■ V p may be 
obtained as an appropriate combination of the vectors Y™{p,x) and X™(p, X ). This completes the proof 
of Proposition 14.31 □ 



The representation (|4.6p and Proposition 14.31 imply, in virtue of Hormander's theory [ITJ[TniliIll2E]> the 
following result: 

Corollary 4.4. The operators L, L\ d t — L, d t — L^ and 2d t — L ® 1 — 1 ® are hypoellipticQ on M 2 ^ 



i+ x K 2 ^ and R+ x R 2 ^ x M 2 ^, respectively. 

In particular, the hypoellipticity of 2dt ~L®1— 1®^ implies that the transition probabilities of the 
dispersion process p(t), 

Pt(p ,dp) = P t (po,p)dp, (4.15) 

have densities (annihilated by 2dt — L <Ei 1 — 1® L>) that are smooth functions of (t,po,p) for t > and 
away from the origin in M. 2d . 



differential operator D on a domain f2 is hypoelliptic if for all distributions /, g such that 
Df = g, smoothness of g on an open subset [/Cfl implies smoothness of / on U. 
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5. CONTROL THEORY AND IRREDUCIBILITY 

The additional important property of the process pit) restricted to R 2 ^ is its irreducibility assured by 
the strict positivity the smooth transition probability densities Pt(po,p) f° r all £ > and po 7^ ^ p. 
The latter property results, according to Stroock-Varadhan's Support Theorem [32], see also (30], from the 
controllability of the process p(t) on R 2 ^ that is established in the following 

Proposition 5.1. For every T > and po ^ 7^ p\ there exists a piecewise smooth curve [0, T] 3 t i-» 
(u™ it)) £ R d2 such that the solution of the ODE 

p = X (p) + J2<(t)Xn(P) (5-1) 

with the initial condition p(0) = po satisfies p(T) = p%. 

Proof. First suppose that po 7^ 7^ p\. Let [0,T] 3 1 H> p(£) be any curve such that 

p(0) = po , rp(0) = Xo , 

p(T) = Pl , rp(T) = xi (5.2) 

and such that p(t) 7^ for all i G [0,T]. Set x(*) = T P{t)- Let 

0(t) = rx(t)+x(t). (5.3) 



Then the formula 



where now 



= 1 7 (ap n p m + W"), (5-4) 



£ + PJ_0 fl= J_ ^ 

+ F + G)F (p 2 ) 2 ' P Fp 2 ' 1 ' 

defines smooth control functions [0, T] 3 t M- (u^(£)) such that Eq. (|5.1|) holds. 
Now suppose that po = 7^ pi. Choose < e < and for < £ < e, set 

p(£) = (l-e-^)xo (5.6) 

and x(t) = rp(t). Then p(£) = (p(£),x(£)) satisfies Eq. with u™(£) = for < t < e, with the 
correct initial condition at t = 0. Note that 

(p(e),x(e)) = ((1 - e-*)xo, e-*xo) ■ (5.7) 

Since, by the assumptions, Xo 7^ 0, we infer that p(e) 7^ and the solution of Eq. (|5.1[) for e <t <T may 
be constructed as in the previous point but taking (|5.7p as the initial conditions at t = e. 

Similarly, if p 7^ = pi then set for T - e < £ < T 



T-t 



pit) = (1-e— )xi, (5.8) 

and x(i) = rp(*)- Then p(£) = (p(£),x(£)) satisfies Eq. (JSTTJ) with u™ (£) = for T — e < t <T , with the 
correct final condition at t = T. One has 

(p(T-e),x(T-e)) = ((1 - e*) X i, e**i) ■ (5-9) 

Since, by the assumptions, Xi 7^ now, we infer that p(T — e) 7^ and the solution of Eq. (|5.1[) for 
< £ < T — e with p(£) 7^ may be constructed as in the first point but taking (|5.9I) as the final condition 
at t = T — e. 

Finally, if po = = p\ , we combine the above solutions for < t < e and T — e < t < T with vanishing 
ujjj with the a solution with pit) 7^ and appropriate u„(£) for e < £ < T — £. □ 

Remark 5.2. Note that the solution p(£) of the ODE (pTT]) satisfying p(0) = p 7^ and p(T) = pi 7^ is 
everywhere nonzero. 
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6. PROJECTION OF THE DISPERSION TO S 2rf 1 
The generator L of the process commutes with the multiplicative action of R+ on R 2d given by 

p ap (6.1) 
for a > 0. It follows that if p(0) 7^ then the projection 

\p(t)} EE 7T(i) (6.2) 

of the process pit) on the quotient space R 2 / /R+ is also a Markov process whose generator may be identified 
with L acting on functions on R 2 ^ that are homogeneous of degree zero. The quotient space R 2 i u /R+ may 
be naturally identified with the sphere 

S 2 "- 1 = {( P ,X) I P 2 +X 2 =R 2 } (6.3) 

for a fixed R and we shall often use this identification below. The transition probabilities Pt(^o', dir) of the 
process 7r(i) are obtained by projecting the original transition probabilities from R 2 ^ to the quotient space. 
Note that the vector fields Xo, X™, Y™, Z™ also commute with the action R + so may be identified with 
vector fields on R 2 / /R+ and Eq. gU) still holds. Viewed as vector fields on S 2 ^ 1 , X™, Y™ and Z™ still 
span at each point the tangent space to S' 2d_1 . It follows that the operators L, L\ dt — L, dt — L\ 2dt — 
L <£> 1 — 1 (g) lJ (with the adjoints defined now with respect to an arbitrary measure with smooth positive 
density on S 2d ~ l , e.g. the normalized standard S f O(2d)-invariant one ^{dir)) are still hypoelliptic and the 
transition probabilities of the projected process have smooth densities P t (7ro;7r) with respect to no(dir) for 
t > 0. Consequently, the process Tz(t) is strongly Feller: for bounded measurable functions / on S 
the functions 

(T t /)(7r ) = J P t (7r ;d7r)/(7r) = j P t (7r ; tt) /(tt) MoW (6.4) 

are continuous (and even smooth) for t > 0. Besides, the projected process is still irreducible since 
P f (7r ;7r) > for all i > and 7r ,7T e 5 2d_1 . The latter property follows from the relation between 
P t (7T ; tt) and PtipolP) an d from the strict positivity of the latter away from the origin of R 2d . 

The gain from projecting the process p(t) to the compact space 5 2d_1 is that the projected process 7r(i) 
has necessarily invariant probability measures /x(d7r). In particular, each weak-topology accumulation point 
for T — ¥ 00 of the Cesaro means 

T 

T- 1 J P t (Tv ;dn)dt (6.5) 


provides such a measure^. Since the (a priori distributional) density n(7r) of an invariant measure is 
annihilated by L\ the hypoellipticity of the latter operator assures that n(7r) is a smooth function. The 
invariance relation 

J Pt{^o,^)n{7v ) fi {dn Q ) = n(n) (6.6) 

S 2d-1 

together with the strict positivity of P t (7To,7r) implies then the strict positivity of the density n(7r) of the 
invariant measure and, in turn, the uniqueness of the latter (different ergodic invariant measures have to have 
disjoint supports, so that there may be only one such measure), see e.g. [30] for more details. One obtains 
this way 

Theorem 6.1. The projected process if(t) has a unique invariant probability measure fi(dw) with a smooth 
strictly positive density n(ir). 

The smoothness of the densities Pt (7To;7r) implies by the Arzela-Ascoli Theorem that the operators of 
the semigroup T t on the space C(S' 2d_1 ) of continuous function on 5 2d ~ 1 with the sup- norm, defined by 
Eq. (16.41) . are compact for t > 0. The uniqueness of the invariant measure implies then that the spectrum 
of T t is strictly inside the unit disk except for the geometrically simple eigenvalue 1 corresponding to the 
constant eigenfunctions, see 30 . It follows that the process 7r(i) is exponentially mixing: 



Probability measures on a compact space form a compact set in weak topology. 
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Theorem 6.2. 



(/i(7r(ti))/ 2 (7r(i 2 ))) t ^t f /iWAi(dTr) f h(n)fi(d-K) (6.7) 
exponentially fast for continuous functions fi , fa . 

7. PROPERTIES OF THE INVARIANT MEASURE 

Due to the isotropy of the covariance (|2.5[) . the generator L of the process 7r(t) commutes with the action 
of the rotation group SO(d) induced on S 2 ^ 1 by the mappings 

(p, X ) & (Op, Ox) (7.1) 

for O € SO(d). As a consequence, the process 7r(i) stays Markov when projected to the quotient space 
Pd = S 2 ^ 1 1 'SO(d). The unique invariant measure fi(dw) of the process n(t) has to be also invariant 
under SO(d) and its projection to Pd provides the unique invariant probability measure of the projected 
proces^. The projected invariant measure may be expressed in terms of invariants of the 5*0((i)-action. 
Such invariants will be chosen as the following dimensionless combinations: 

• for d — 1 where Pi = S 1 

* = - , (7-2) 
P 



• for d = 2 where P 2 = PC 1 , 



px , p 1 x 2 -p 2 x 1 ,~„\ 

— — and y = - (7.3) 

P P 



with z = x + iy providing the inhomogeneous complex coordinate of PC 1 , 
• for d > 3, 



p x A Vp 2 x 2 - (p-x) 2 

x = — 5- and y = ^ . (7.4) 

p- p l 

Note that the right hand side of the d > 3 expression for y would give in d — 2 the absolute value of y. 
The quotient spaces Pd are not smooth for d > 3. 



7.1. d — 1 case. In one dimension, Eq. (12.8[) implies that 

dx = -±(x + x 2 )dt + dS{t). (7.5) 

The invariant probability measure on S 1 is easily found [141133] to have the form dfi — n(x)dx with 

X 

n(x) = Z- 1 (^-^T2Bj(i x3+x2 ) J c^T^^ x ' 3+x ' 2 )dx^dx, (7.6) 

— oo 

where Z is the normalization constant. Since the normalized rotationally invariant measure on S 1 = 
{(p, x) I P 2 + X 2 = R 2 } has the form d[i Q = -^jf^, it follows from our general result that the density 
n(x) = 7r(l + x 2 )n(x) of the invariant measure relative to dpo must be smooth and positive at x — oo, i.e. 
at the origin when expressed in the variable x^ 1 . In particular, 

r)(x) = 0(\x\- 2 ) for |x|->oo, (7.7) 

which may also be easily checked directly. 

In one dimension, the generator L given by Eq. (I2.1ip acts on a function f(x) according to the formula: 

(Lf)(x) = -l(x 2 + x)d x f{x) + \{A + 2B)dlf{x). (7.8) 



It coincides with the generator of the process satisfying the SDE (|7.5[) . The trajectories of the latter process 
with probability one explode to — oo in finite time but, in the version of the process that describes the 
projectivized dispersion of the one-dimensional inertial particle, they re-enter immediately from +oo. 



^To see the uniqueness, note that averaging over the action of SO(d) maps C(S 2d 1 ) to C(Pd) and that dual map sends 
invariant measures for the projected process to invariant measures of n(t). 
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7.2. d = 2 case. 

In two dimensions, the invariant measure on S 3 has to have the form 

d V = ^v( z ,z)d 2 zdaTg(p) . (7.9) 
On the other hand, the SO (4)-invariant normalized measure on S 3 is 

= J^Vo(z,z)d 2 zdarg(p) (7.10) 



Vo{z,z) = ; ,_ |2 . 2 ■ (7.11) 



with 

1 

7T(1 ~|Z| 2 ) 

It follows from the general result obtained above that the density of dfi relative to d/io 

n(z,z) = (7.12) 

has to extend to a smooth positive function on PC 1 , i.e. to be smooth and positive at zero when expressed 
in the variables (z ,z ). In particular, 

T](z,z) = 0{\z\- 4 ) for |*|->oo. (7.13) 

The unique invariant probability measure of the Markov process obtained by projecting 7r(i) from S 3 to 
S 3 /SO(2) = PC 1 has the form (|7.9p with ■^da,vg(p) on the right hand side dropped. Note that the relation 
(|7.13p implies that 

oo 

r}(x,y)dy = 0(M~ 3 ) for \x\ ->• oo . (7.14) 



by changing variables y t— > v 1 + x 2 y in the integral. Such behavior was heuristically argued for and 
numerically checked in [2]. 

In two dimensions, the generator L of Eg. (|2.11l) acts on SO(d) invariant functions f(x,y) according to 
the formula: 

(Lf)(x,y) = -±(x 2 -y 2 + x)d x f(x,y) 
--(2xy + y)d y f{x,y) 

+ \(A + 2B)d 2 J{x,y) + \Ad 2 y f{x,y). (7.15) 

It coincides with with the generator of the process z(t) = (x + iy)(t) in the complex plane given by the 
SDE M 



dz = -^(z + z 2 )dt + VA + 2Bdf3 1 (t) + iVAd/3 2 (t), (7.16) 
where /5 ' (f) and /3 2 (t) are two independent Brownian motions. 
7.3. d > 3 case. 

Finally, in three or more dimensions, the invariant measure on S' 2d_1 has to have the form 

dfi = r)(x,y) dxdy d[0] , (7-17) 

where O 6 SO(d) is the rotation matrix such that 0~ x p is along the first positive half- axis in R d and 
0~ 1 x lies in the half-plane spanned by the first axis and the second positive half-axis. Note that, generically, 
O is determined modulo rotations in (d — 2) remaining directions. d[0] stands for the normalized SO(d)- 
invariant measure on SO(d) / SO(d — 2). In the same notation, the 5 , 0(2<i)-invariant normalized measure on 
S 2 ^ 1 takes the form 

d[io = Vo( x i y) dxdy d[0] . (7-18) 

for 

(d-l)2 d 'V- 2 

Vo{x,y) = —7— — 2 2 . d , 7.19 
7r(l + x z + y l ) a 

as is shown in Appendix |Aj As before, it follows from the general analysis that the function 

/ v V(x,y) 

n(x,y) = —, r (7.20 

Vo{x,y) 
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is smooth and positive on the sphere S 2d 1 = {(p, x) I P 2 + X 2 = R 2 }- In particular, this implies that 

r)(x,y) = 0(y d - 2 ) for y\0 (7.21) 

i.e. for p and x becoming parallel or x 2 becoming small and 

T](x,y) = 0{\x\- 2d ) for \x\ -> oo (7.22) 

when p 2 — > but the angle between p and x stays away from a multiple of ^. The smoothness and 
positivity of n(x,y) on S 2 ^ 1 imply (again by changing variables y H> y/l + x 2 y in the integral) that now 

rj(x,y)dy = G{\x\' d ^) for |x|->oo. (7.23) 



o 

A straightforward calculation shows that, in three or more dimensions, the action on L on S f O(<i)-invariant 
functions f(x, y) is given by a generalization of Eq. (|7.15[) : 

(Lf)(x,y) = -L(x 2 -y* + x )d x f{x,y) 

-L(2xy + y -^V)d y f(x,y) 

+ l(A + 2B)d 2 J(x,y) + \Ad 2 y f{x,y). (7.24) 
It coincides with the generator of the process z(t) = (x + iy)(t) in the complex plane given by the SDE [2|IT7] 
dz = -~(z + z 2 -i l£±zQ. )dt + VA + 2Bdf3 1 (t) + iVJdf3 2 (t) (7.25) 
which upon setting d = 2 reduces to the SDE (|7.16[) . 

8. ABSENCE OF EXPLOSION IN THE COMPLEX (HALF-)PLANE 

Let us set 

Qd -\m + if d>3, (8 - 1} 

where 

m+ = {(x,y) \y>0} 

is the upper-half plane. Note that Qd may be identified with an open dense subset of the quotient space 
P d = S 2d ~ 1 /SO(d) using the 5'0(d)-invariants (|7T3f or (jTH) on S 24 " 1 . We shall often use the complex 
combination x + iy as a coordinate on Q^- 

In the present section, we shall show that for d > 2 the unique solution of the SDE (17.251) starting from 
z G Qd remains in Qd for all times t > with probability one. This will also have to be the property of 
the projection of the process n(t) to the quotient space Pd = 5 ,2d_1 / SO(d) when described in the complex 
coordinate z = x + iy. Indeed, the coincidence of the generators of the two processes will assure that they 
have the same law. Let us start by generalizing and simplifying (|7.25|) . 

Let w(t) — z(t) + 1/2, where z{t) solves (|7.25l) with z(0) = x + iy £ Qd- Clearly, w(t) satisfies an SDE 
of the form 

dw = i (-w 2 + a + i^g!) dt + J^dP^t) + t] ^dl3 2 (t), (8.2) 

where a = ai + ia 2 G C, b > 0, K\ > 0, k 2 > 0, and ^(t) and (3 2 {t) are two independent Brownian motions. 
When d = 2, the term proportional to b(d — 2) is absent from (I8.2[) . When d > 3, we suppose that 
rb(d— 2) > K2- Clearly in (17.25[) . all of these assumptions are met under the given substitution. Since w(t) 
is a horizontal shift of z(t), w{t) stays in Qd with probability one for all times if and only if z(t) does. 
Employing methods of refs. [15 1 I26 1 I30] . we shall estimate the time at which the process w(t) leaves Qd- To 
this end, it is easy to see that there exists a sequence of precompact open subsets {O n \ n G N} of Qd such 
that O n t Qd as n — > oo. Thus we may define stopping times: 

r„ = inf{s > | w(s) G O c n }, (8.3) 

for n G N. Let be the finite or infinite limit of r„ as n — > oo. 

Definition 8.1. We say that the solution w(t) is non-explosive if 

P[r 00 = oo] = l. (8.4) 
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Naturally, in order to show that w(t) remains in Qd for t > with probability one, it is enough to prove 
that w(t) is non-explosive. 

Let M be the generator of the process w(t) = x(t) + iyt- We see that for / G C°°(Qd) ■ 

(M d f)(x,y) = - 1 7 (x 2 -ij 2 -a 1 )d x f(x,y)~^(2xy-a 2 -rb(d-2)y- 1 )d y f(x 1 y) (8.5) 

+ ^d*f(x,y) + ^d 2 y f(x,y), 

where the term rb(d — 2)y _1 is absent for d — 2. Let us define 

oo if d = 2 

{(x, y) G M 2 | y = 0} U oo if d > 3 



d Qd={<, , ^ _ 2 , !P L„' (8-6) 



with oo denoting the point compactifying R 2 . To ensure condition (I8.4p . it suffices to construct a (Lyapunov) 
function $d € C°°(Qd) that satisfies: 
(I) ®d(x,y) > for all (x,y) G , 
(II) $ d (x,y)->'Oo as (x, y) -> 9Q d , (x,y)eQ d , 

(III) Ald&d(x,y) < C$d(x,y) for all (a;,y) G Qd, where C > is a positive constant. 
See, for example, [2B]. We will show: 

Theorem 8.2. If n x > and rb(d - 2) > k 2 > then there exists $ d G C°°(Qd) iftat satisfies jl|, dTTJ) . 
and 

M d $ d (a;, y) -> -oo as (x, y) -> <9Qd, (a;, y) G Qd- (IV) 
Given such <E> d , clearly $d + 1 will satisfy (jU, dTTJ) and pIT|) . We will then have: 

Theorem 8.3. Under the assumptions of Theorem \8. 6 A the solution w(t) of the SDE (|8.2I) stays in Qd for 
all times t > zwi/i probability one if w(0) = x + iy G Qd- 

Corollary 8.4. This implies the same result about the solution z(t) of the SDE (|7.25[) with A > and 
A + 2B > 0. 

The existence of the Lyapunov function with the properties asserted in Theorem 18.21 has another conse- 
quence. It allows to show that 

lim liminf i / P t {w,O c n )dt = (8.7) 

n— >oo T— yoo 1 Jq 

for the SDE (|8.2I) and O n f Qd as before, implying the existence of an invariant measure on Q d , see Theorems 
4.1 and 5.1 in Chapter III of [15]. If the generator of the process is elliptic, then the same tools that we used 
for the projectivized dispersion (i.e. hypoellipticity and control theory |30| ) show that the invariant measure 
must have a smooth strictly positive density and be unique. This gives: 

Theorem 8.5. Under the assumptions of Theorem ] 8. l A the system (|8.2[) on Qd has an invariant measure 
which is unique and has a smooth strictly positive density if K\ > 0. 

Remark 8.6. Theorem 18.51 allows to reaffirm and strengthen what has already been proven earlier since 
it implies the existence of an invariant measure for the system (I7.25|) if A > and A + 2B > and its 
uniqueness if A + 2B > 0. Given the non-explosivity result of Corollary 18.41 the approach taken earlier 
implied the existence and the uniqueness of an invariant measure for the system (|7.25p under more stringent 
conditions: A > 0, A > \B\ and A + (d + l)B > 0. 

The construction of the Lyapunov function $d is split up into two cases: d = 2 and d > 3. The existence 
of <&d for d > 3 will be easy, given <J> 2 . Thus we shall first construct $ 2 - 

8.1. d = 2 case. 

It is not easy to write down a globally defined function $2 that satisfies (P), ©, an d (|IV|) in all of Q2 = K 2 . 
This is because the signs of the coefficients of the vector fields in M 2 vary over different regions in R 2 . We 
shall thus construct functions that satisfy these properties in different regions, the union of which is M 2 . We 
shall then glue together these functions to form one single globally defined $2- One should note that this 
idea is similar in spirit to that of M. Scheutzow in [21]. Let r — \J 1 x 2 + y 2 . For the rest of Subsection 18.11 
we will drop the use of the subscript 2 in M 2 and $2- We first need the following: 

Definition 8.7. Let X C K 2 be unbounded. We say that a function f(x,y) — > ±00 as r — > 00 in X if 

f(x, y) -> ±00 as (x, y) -> 00, (x, y) G X. 
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Definition 8.8. Let Id 2 be unbounded and let (p E C°°{X) satisfy 

(i) f > for all (x, y) € X, 

(ii) (p — > oo as r — > oo in X, 

(iii) M(p — > — oo as r — > oo in X. 

We call </? a Lyapunov function in X corresponding to M and denote 

Af(a, Ki, k 2 , X) — { Lyapunov functions in X corresponding to M } . 
We shall abbreviate "Lyapunov function" by LF. 

Definition 8.9. Let X C R 2 be unbounded and /, g : X — > R. We shall say that / is asymptotically 
equivalent to g in X and write / ~x g if 

f fx, w) 
lim J -j-^- = 1 , 

where the limit is taken only over points (x, y) € X. 

It is clearly sufficient to construct LFs in regions that cover R 2 , except, possibly, a large ball about 
the origin. The constructions will be done in a series of propositions. The possibly daunting multitude of 
parameters is designed to make the gluing possible. There is a total of five LFs in five different regions and 
the details that follow are not difficult to verify. The crucial LF is the fifth one, ip^, defined in a region where 
explosion occurs in a nonrandom equation, i.e., when a = k± = k 2 = in (|8.2p . 

Proposition 8.10. Let X 1 = {x > 1} C R 2 , C x > 0, and 5 E (0, 1/2). Define 

Mx,y) = C 1 (x 2 +y 2 ) s / 4 . (LF1) 

We claim that ip% g J\f(a, ki,K2, Xi) for all a £ C, Ki > 0, k 2 > 0. 

Proof, ipi is nonnegative everywhere in R 2 , hence everywhere in X±. ip — > oo as r — > oo in all of R 2 , 
hence in all of X%. It is easy to check that d xx p>\ and d yy <pi both go to zero as r — > oo. Thus dropping 
second order terms in the expression for Mipi, we have 

tM ^ -^x(x 2 + » 2 )^ 4 + ^ ^ + *2y 

^ Xl _^a;(x 2 +y 2 )^ ^ -oo (8.8) 
as r — > oo in X\, since x > 1 in X%. □ 

We need a remark before we move onto the next region. Let R C R 2 be the real axis. 
Remark 8.11. Let f(x,y) — u(x, \y\) be a twice differentiable function in X \ R. Then 
(rM/)(x,t/) = Kiu xx (x,\y\) + K2U\y\\ y \(x,\y\) 

+ (y 2 - x 2 + ax)u x (x, \y\) + {-2x\y\ + sgn(y)a 2 ) M| y |(x, \y\). 
Proof. Apply the chain rule to the operator d\ y \. □ 

In the following arguments, often the function will be of the form f(x, y) = it(x, \y\). The above remark will 
allow for simplifications in the argument for property (iii) in Definition 18.81 

Proposition 8.12. Let C 2 > 0, 6 E (0, 1/2) and 

V2 = C 2 (-x+\y\ s / 2 ). (LF2) 

Then ip2 & N{a, K\, K2, X2) for all a E C and all k\ > 0, K2 > 0, where 

X 2 = {-2 < x < 2} n {|y| > 2 2 ' 5 }. 

Proof. tf2 is indeed smooth in X2 since X2 is bounded away from R. Note that the region was chosen so 

that if2 > in Moreover, since x is bounded in this region, r — > 00 in X2 if and only if \y\ — > 00. 

Hence, if2 — > 00 in X2- By Remark 18.111 and noting that d xx Lp 2 — and that d\ y \\ y \Lp2 — > as |y| — > 00, 
we have 

rM^ 2 (x,y) ~x 2 C 2 (x 2 - 2 ; 2 -a 1 ) + C 2 -(-2x| 2 /|+sgn( 2 /)a 2 )|y|^ 1 
— x 2 -C 2 y 2 -> -00 

as r — > 00 in X 2 . □ 
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Proposition 8.13. Let C 3 > and <5 e (0,1/2). Define 

^ = C 3 (?±g- ) (LF3) 

on X 3 = {x < -1} n {\y\ > 1}. Then ip 3 e A/"(a, ki,k 2 , -^3) for all a e C and all m > 0, k 2 > 0. 

Proof. Smoothness of ip 3 is not a problem in this region as we are bounded away from R in X3. Clearly, 
¥3 > and note that ^-foo as r — > 00 in X 3 . After dropping the 5(5 — l)-terms which are negative, we 
obtain: 



TM<p 3 ( Xl y) < C 3 6 



x" + y \ 2 / 15x 2 1 \ x 3 



|y|3/2 



«1 , ,,,, + K2 ,- /Q - Ti iTTo + 



y |3/2 \^4|y| 7 /2 4|y| 3 / 2 / |y| 3 / 2 



+ %| 7 + f-p^ + sgn(y)a 2 



y |3/2 *» - ^ 2 | y | 5/2 2 | y | 1/2 

= <5xy>3 — > —00 (8.9) 
as r — > 00 in X 3 since x < — 1 in X 3 . □ 

Proposition 8.14. Let C 4 > 0, 77 > 1 and <5 e (0, 1/2). Define 

I 125 , I |2<5 

¥>4(s,y) = C4 FI , Vr (LF4) 



on 



A 4 = {x < -1} n I ry^f^+l)-^ < |y| < 2 J . 
Then 934 € A/"(a, /ci, k 2 , Xi) for all aeC and all /ti > 0, k 2 > 0. 

Proof. Note that <p 4 is smooth in X4 since this region excludes both x and y axes. Moreover, 954 — > 00 as 
r — »■ 00 in I4 since then x must approach 00 as r —¥ 00, and ?/ is bounded above. Dropping insignificant 
terms in the expression for Mip^, we see that in X 4 : 

. I r |2<5+1 „ , N | r |2<5 

rMipi(x, y) < C 4 (-5^ + K2 ls(l5 + l)±^+6\x\\y\^ 
~ ai2( * , ,3, + sgn(y)a 2 ^ 2 - sgn(y)a 2 -<5- 



|y|2» - v- / |y|3*+2 



/ . . 3/3 

Ivl 



I 125+1 

< -CW^-l/^ 2 )'-' -> -00 

in X4 as r — > 00. □ 

Proposition 8.15. Let C 5 , (3 > and £ > such that 2k 2 > E(3, let £ > 1, and let 

<p 5 (x,y) = C 5 (E\x\P-y 2 \x\f } + 1 ) (LF5) 

be defined on 

* 5 - {*<-i}n ||y|<iyZ|. 

Then 935 € A/"(a, «i, k 2 , X 5 ) for all aeC and all Ki > 0. 
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Proof. The fact that ip5 is smooth in X§ is clear as x < — 1 in X§. Again by the choice of Xs, ^5 > 
and <p$ — > 00 as r — > 00 in X§. Dropping irrelevant terms in the expression for Mips, we see that in A5: 

rM<p 5 (x,y) < C 5 (K 1 E/3(P-l)\xf- 2 -2K 2 \xf +1 +Ep\xf +1 + (f3 + l)y 4 \xf) 

+ C 5 {a^-Eftxf- 1 + ((3 + l)y 2 \xf) - 2a 2 y\xf +1 ) 

~x 5 C 5 (Ep - 2 K2 ) \x\ fj+1 -> -cm 

as r — > 00 in X5, as |i| must approach 00 when r — > 00 in X5. □ 

We now have our desired LFs. It is not obvious, however, that the regions Xi, X 2 , . . . ,X& cover R 2 
except, possibly, a bounded region about the origin. To assure that one has to show that X4 and X§ 
overlap. In order to make this more tangible, we will choose some of the parameters given in the previous 
propositions. With the choices that follow, however, we first need a lemma that says that varying the diffusion 
coefficients (ki, k 2 ) is permitted. This lemma will also be of crucial use later when we glue the LFs to form 
a globally defined 

Lemma 8.16. Fix k 2 > and suppose that $ G Af(a, k%, K2,R 2 ) for all aeC and all k% > 0. Then for 
every t 2 > 0, we can find \1/ G Af(a, i\, /<2,R 2 ) for all a G C and all i\ > 0. 



For the proof of Lemma [8T6l we temporarily use the notation M? Ki R \ for the generator M given by (|8.5|) . 

Proof. Let 77 > be such that rf t 2 = n 2 . Define ^(x,y) — ^(rjx,rjy). For the function "J, smoothness 
and properties (|T| and (JTTJ) are immediate. Let s — rjx and t = r\y. Then by the chain rule 

tM? , *V(x,y) = -M? v \ 3 Ms,t) 

= -M^\ v«(s,i) -> -00, 



as r = y/x 2 +y 2 — > 00 . □ 

By Lemma T8. 161 it is enough to find a function $ 6 Af(a, K\, k 2 , R 2 ) for some fixed K2 > for all a G C 
and all ki > 0. All of the <fi satisfy these criteria. In fact, tpi for i — 1,2,3,4 work more generally. The 
reason that cps only works for 2^2 > Ej3 is due to the fact that when a = K\ = n 2 — 0, the solution to 
(|8.2p has an explosive trajectory along the negative real axis. 

Now we choose some parameters. Let E = 5 and £ = \/5/2, so as to make 

i 5 = {i<-m {MM 1/2 < 2}. 

Let /3 = ^5, 5 = jK 2 , and k 2 G (0, 1). Then Ef3 = §§k 2 < 2k 2 and <5 G (0, ±). Hence for alii = 1, 2, . . . , 5, 
tfi is a LF in the region JQ. Decrease K2 > so that 

X 4 D {x < -1} n {1 < |y||a;| 1/2 < 2}. 

One can easily check that X4 and A5 overlap in such a way that we have covered all of R 2 with X\ , X2 , . . . , X$ 
except a bounded region about the origin. 

We fix k 2 > sufficiently small (this will be made precise later). We will construct a function <E> G 
Af(a, «i, K2,R 2 ) for all a G C and all «i > 0. The idea is as follows. Note that ipi and ip 2 are LFs in 
the region X\ n A 2 . We shall define a nonnegative smooth auxiliary function £(x) G C°°(R) such that 
C(x) =0 for x > 2 and £(x) = 1 for a; < 1, satisfying some additional properties. We will then show that 
the combination 

(1 - ()ipi + c,ip 2 

is a LF in the larger region X\ U A2. Proceeding inductively this way, we shall construct a LF in all of M 2 . 

Let us first define some auxiliary functions needed to construct such a $. Let ( : R — > R + be a C°° 
function such that 



CO*) 



for x < 1 , 
for x > 2 , 



and C'C 2 -) < for all x G (1,2). We define the smooth function fj, : R — > K + as the horizontal shift of £, 
three units to the left, i.e., 

fi{x) = C(x + 3) for iGl. 
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Let 

'C(M) for a; < -2, 
v(x,y) = I for \y\ > 2, 

for x > — 1 , 

and assume that v is C°° outside of the ball B4. Let q : (—00, —1] xR- 

1 if |x| 1 /2|y|> 2 , 



be defined by 



and 



Let 



q{x,y) = 



r(f) 



Sl^yl -1 if 1< ^l^lyl < 2 , 

if \ x n\<i, 



exp 




( l-(2t-l) 



if < < < 1 
otherwise . 



s (x) = i-/ r{t)dt. 
^ J —00 

where AT = J R r(t)dt. Now define a function on R 2 by 



p(x,y) 



s(q(x,y)) if ap < — l, 
1 if M>3, 

1 if x>-l/2. 



Clearly, p is C°°(R 2 ) outside of B 4 



Let r(5) = max + 4^. Define tp for a; 2 + y 2 > r 2 (5) by 



<P1 






if 


x> 2, 


(<P2 H 


-(1 


-CM 


if 


1 < x < 2, 


¥>2 






if 


- 1 < x < 1, 




Ki 


- 


if 


- 2 < x < -1 , 


</>3 






if 


x < -2, |y| > 2 


1/(^4 - 




- ^)¥>3 


if 


x < -2, 1 < \y\ < 2, 


</>4 






if 


x<-2, |]/| < 1, \x\ 1 / 2 \y\>2, 


W>4 H 


-(i 


- P)^5 


if 


x<-2, 1< Ixl 1 / 2 ^! <2, 


¥>5 






if 


x < -2, Ixl 1 / 2 ^! < 1, 



and now 



$(x,y) 



if 



+ 2/ 2 > B r{s) 



1 arbitrary positive and smooth if x 2 + y 2 < B r /$\ . 
It is easy to see $ can be chosen to be nonnegative and C°°(R 2 ). With the aid of Lemma [8.161 the following 



lemma implies Theorem 18.21 in the d — 2 case. 

Lemma 8.17. For K2 sufficiently small, $ £ Af(a, K\, M 2 ) for all a £ C and all K\ > 0. 

Proof. Clearly, $ is smooth and satisfies properties I and II. Since M(fi — > —00 as r — > 00 in Xi for each 
i, all we must verify is that M$ — > — 00 as r — > 00 in the overlapping regions. Let us recall the choices that 
have already been made: 



E 



2 ' r 4 

and note that n 2 £ (0,1) was chosen such that 

Xa D {x < -1} n {1 < |y|M 1/2 < 2} 
Pick C 5 > Ci — C 3 > C 2 > G\. Consider first 

ipi ■= C<^2 + (1 - CM 

defined in the region 

Yi = {Kx<2}n B c r{S) . 
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We have 

TM-01 = (TMV2 + (1 - tMi P1 + (y 2 - X 2 + Cli)C(tp 2 - tfil) 
+ Kl(£'((p2 - Ifil) + 2('(d x (p 2 ~ d x ipi)) 

-(C 2 y 2 -(l-0^x(x 2 +y 2 ) s / 4 + (y 2 -x 2 + a 1 )C(C 2 -C 1 )\y\ s / 2 (8.10) 

+ K 1 (C'(C 2 -C 1 )\y\ s / 2 -2('C 2 ). 

Note that if x is bounded away from 1 and 2 in (1,2), the dominant term above is y 2 ('(x)(C 2 — Ci)\y\ 5 / 2 — > 
— oo. Note also that as x — > 1 or x — > 2, £" — > 0. ,But, the first two terms in the expression above decay 
to — oo at least as fast as — C\y\ 5 / 2 , where C > is a constant independent of £ and x. Thus we may 
choose e > so that whenever x € (1, 2) \ (1 + e, 2 — e), Mipi decays at least as fast as — D\y\ 6 / 2 , where D 
is some positive constant. 
We now consider 

ip2 ■= + (1 - M)^2 

in the region 

Y 2 = {-2<x<-l} n B c r(s) . 

We have 

tM^ 2 = [iTMip 3 + (l- n)rMip 2 + (y 2 -x 2 +ai)n'(ip 3 - tp 2 ) 

+ Ki(//'(¥> 3 - ip 2 ) + 2// (d x tp 3 - d x ip 2 ) 
~y 2 »8x<p 3 - (1 - »)C 2 y 2 + (y 2 - x 2 + ai)// (C 3 - C 2 )\y\ s/2 

+ K 1 (p"(C 3 -C 2 )\y\ s / 2 +2fi'C 2 ) 
^ Y2 n C 3 x\y\ s / 2 - (1 - n)C 2 y 2 + (y 2 - x 2 + ai )^(C 3 - C 2 )\y\ s ' 2 

+ K 1 (p"(C 3 -C 2 )\y\ s / 2 +2fi'C 2 ). 

Note that, for the very same reasons as in the case of ip\, Mip 2 — > — oo as r — > oo in F 2 - 
Let 

V>3 := v*P4 + (1 - f)<£3 

in the region 

y = {x<-2} n {i< | V | <2} n s r c (5) . 

Thus 

rM-03 = vrMtpi + (1 - v)tMlp 3 + (-2x\y\ + $,gn(y)a 2 )(d\ v \v)(tp i - ip 3 ) 
+ K 2 (d 2 yl v)(Lp 4 - ip 3 ) + 2K 2 (d\ v \u)d\ y \(<fi4 - Vz) 
~y 3 vtMlp4 + (1 — v)rM(p 3 — > — oo . 

This is true since we chose C3 = C4. Hence the first order term 

{-2x\y\ + sgn(y)a 2 )(d\ y \is)(ip4 - <#$) 

approaches infinity at worst as fast as C\x\, where C is a positive constant. The second order term 
K 2{y\ y \v){fPi — <p 3 ) is a bounded function in this region. Moreover the term 2K 2 (d\ y \i>)d\ y \ (934 — ip 3 ) at 

worst approaches infinity as fast as D|x| 25 , where D is some positive constant. But, both Mip 3 and M<p4 
approach negative infinity at least as fast as —D\x\ 2S+1 , where D is another positive constant. This gives 
the desired result. 
Let 

tp 4 ■= m + (1 - p)<P5, 

in the region 

Y 4 = {1< \x\^ 2 \y\ <2}n{x< -2} n B c r(5) . 

Note that 

rMV-4 = s(q)TM<p 4 + (1 - s(q))rMip 5 + d x (s(q))(<p 4 - f 5 )(y 2 - x 2 + ai) 
+ 9\ y \(s(q))((p4 - <p 5 )(-2x\y\ + sgn(y)a 2 ) + K 1 d 2 x {s{q)){^ i - ip 5 ) 
+ 2K 1 d x (s(q))d x (ip 4: - <p 5 ) + K 2 d 2 yl (s(q))(ifi4 - ip 5 ) 
+ 2K 2 d lyl (s(q))d lyl (ip 4 - ip 5 ) . 
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Note that in the expression above, we may drop the Kid^(s(q))(ip4 — ips) and 2Kid x (s(q))d x (ip4 — ^5) terms, 
as they are asymptotically less than other terms. Dropping other obviously insignificant terms, we obtain: 

tMt/> 4 ^y 4 s(q)TMip 4: + (I - s(q))rM(pr - d x (s(q))((p4 - ip 5 )x 2 
- 9\y\(s(q))(<f4 - (fi5)2x\y\ + K 2 d 2 yl (s(q))(ip 4 - </? 5 ) 
+ 2n 2 d\y\{s{q))d\ v \ (ip 4 - tp 5 ) 
< F(x,y), 
where F is a smooth function satisfying: 

F ~ n -s( q )C A (l-l/ V 2 )^\x\^ s+1 -(l-s( q ))^\x\^ s + 1 (8.11) 
- d x (s(q))(ip4 - (p 5 )x 2 - d lyl {s(q))((p 4 - f 5 )2x\y\ 
+ K 2df y \(s(q))(ip 4 - ¥>s) + 2K2d\y\(s(q))d\ y \(ip 4 - (p 5 ) . 

By the choice of C5 > C4, it is easy to see that for large \x\ in this region, there are constants C, D > 
such that 

III 

ted\y\(s(q))(<P4 ~ <P5) + 2K 2 d\ y \{s(q))d\ y \{(p4 - ip 5 ) < Cn 2 r{q) 

y\ x i y) 

and 

-d x (s{q))(ip 4 - ^ 5 )x 2 - d lyl (s{q))( m ~ V5 )2x\y\ < -Dr(q)\x\^ 5+1 , 

where the constants C, D and the function r(q) are independent of k 2 . r(q) is a function that goes to zero 
faster than any power of the function 

g(x,y) :H(M 1/2 M-1)(M 1/2 M-2)| 
as Ixl 1 / 2 !?/! — ¥ 1 or 2. But note that, for all e > 0, we may choose k 2 so small so that 

Ck 2 



D > 



g(x,y) 2 



for 1 + e < I 1 1 / 2 1 2/ 1 < 2 — e. From the first two terms in (|8 . 11 1) . we obtain at least — C'S x^ s+1 decay for 
some C > independent of k 2 for all 1 < x 1,/2 |y| < 2. But since every other term in the expression goes 
to zero faster than every power of g as | rzr | 1 / 2 1 2/ 1 — ^ 1 or 2, we can choose an e > so small as above so 
that M1P4 — > —00 for all 1 < M 1 / 2 !?/! < 2. This completes the proof. □ 

8.2. d > 3 case. 

Here we shall complete the proof of Theorems 18.21 for d > 3. Recall that for d = 2 and e > sufficiently 
small, we defined $ 2 (x,y) := $(x,y) € C°°(Q 2 ) that satisfies flU, (HO), and (|TV| for all a £ C, k x > 0, and 
k 2 € (0, e). Lemma 18.161 implied then that for k 2 > arbitrary, the function & 2 ^(x,y) := <5> 2 (r]x,r]y) € 
C°°{Q 2 ) satisfied ©, JU), and (JlVj) for all a, provided that e^" 3 = 2n 2 . Let is fix k 2 > 0. For d > 3, we 
shall define 

$d, n ■= $2,,, + log(l + log 2 ( W /2)) . (8.12) 

Lemma 8.18. For fixed k 2 > 0, is a smooth function on Qd that satisfies (JTJ) , (p}, and (JTVj) . 

Proof. By definition, is smooth and nonnegative in Qd = H+. Clearly, Qd.r/ — > 00 as (x,y) — ¥ dQd = 

{(x, y) € R 2 : y = 0} U {00}, (x, y) € Qd- Thus we must verify property (|IV[) . To this end, note that: 

TM d <S>d, v = TM d $ 2 ,n + TM d (log(l + log 2 (to/2))) 

= tM 2 $ 2 , v + r h -^-^d y ^ v + (~2xy 2 + a 2 y + rb(d - 2) - n 2 ) 21 o§(W 2 ) 

V y 2 (I + log 2 {riy/ 2)) 

2^ 2 (l-log 2 ( W /2)) 

2/ 2 (l + log 2 ( TO /2))2' 

Case 1. Suppose first that y > 2r\~ x . It is easy to check that there exist positive constants K\, K 2 > 
such that 

, < K U (a 2 y + ro(d - 2) - n 2 ) + ^-^2)) < 

2/ v2 '"- { 2U y ' 2 % 2 (i + log 2 ( w /2)) y 2 (i + log 2 (vy/2)) 2 - 
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2x - *T < ^ (8.13) 



Ax- — 2 7777 S — -K4X. (°T4) 



If also x > — 2?7 1 then 

l + log 2 ( TO /2) 
for a positive constant K3, whereas for x < —2r/^ 1 , 

1 + log 2 (to/2) 

Since M 2 4>2,? ) — > —00 as (x, y) — > dQ d and y > 2rf x , and, besides, if x < — 2r\~ 1 and y > 2rj~ 1 then 
$2,r){x,y) is equal to the rescaled function ip 3 so that, by (18.91) . M 2 $ 2 ,r) < ifsx(x 2 +J/ 2 ) 15 / 4 for some > 0, 
we infer that M d <$> d>v — > —00 as (x, y) — > dQ d and y > 2vT x . 

Case 2. Suppose now that < y < 2rj~ 1 . If |x| < 2rj~ 1 , then (x, y) — > <9<2d if and only if y | 0. Since 
$2,77 is smooth on K 2 , there exists a constant Kq > such that 

Tb(d - 2) <9<J> 2 n „ 1 

y oy y 

for (x,y) G (— 2t] _1 , 2r)~ 1 ) x (O^^" 1 ). Hence, recalling the assumption rb(d — 2) > K2, we have on this 
rectangle: 

1 , i „ 2 , , ,^ ^ \ 21og( W /2) , 2^(1 - log 2 ( W /2)) 



rM d <f> d>v < K 6 ^ + (-2xy 2 + a 2 y + rb(d - 2) ~ k 2 ) 

< K7 



y V ' y 2 (l + log 2 ( W /2)) y 2 (l + log 2 ( W /2))^ 

1 , 2 K2 (l-log 2 (TO/2)) , 



y y 2 (i + log 2 ( TO /2)) 2 
as y \. 0. 

If x > 277" 1 , then $ 2 ,r; is equal to the rescaled function tpx- We see that by (|8.8[) there exist constants 
-Kg, -Kg i ^10) -^li > such that 

rM d # d ,, < ~K 8 x(x 2 +y 2 ) 5 / 4 +X 9 (x 2 + 2 / 2 ) 5 / 4 - 1 

2Iog( W /2) , 2 K2 (l-log 2 (w/2)) 



+ (-2xy 2 + a 2 y + rb(d - 2) - n 2 ) 



y 2 (l + \og\ m /2)) y 2 (l + log 2 ( W /2)) 2 



y y 2 (l + log^y/2)) 2 

Note that as (x, y) — > SQd in this region then x —¥ 00 or y 4- 0- It is thus easy to see that rM^d.ri —00. 
If x < — 2t7~ 1 , then it is easy to check that d y & 2tV is bounded above by the choice of C3 = C4 and C5 > C4. 
Then, for some i^ia > 0, 

rM d S dv < rM 2 * 2v+ K 12 l + 2 ^ 1 - 1 °(^ 2 » 

,n - ,v y y 2 (i + io g 2 ( w /2)) 2 

so that that TM d $> di „ — > —00 as (x, y) — » in this region. □ 

9. TOP LYAPUNOV EXPONENT 

The Lyapunov exponent A for the dispersion process pit) — (p(t),x(t)) is the asymptotic rate of growth 
in time of the logarithm of the length y/ p 2 + x 2 . Suppose that the process starts at t = from po = 
(POiXo) 7^ 0. Anticipating the existence of the limit below, we shall define: 



lim i((lnVp 2 (T) + X 2 (T) - \n^/p 2 +x 2 )) / 

T 

T lim i/|(lnV?W+^)* 


T 

lim^ 1 J (l In vVW + X 2 ^)) * 


T 

Hm ^ / ( /(£ In VP 2 + X 2 ) Pt(po, dp)) * , (9-1) 
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where L is the generator of the process pit) given by Eq. p. lip . Note that the function L In \J p 1 + x 2 
is smooth on R 2d \ {0} and homogeneous of degree zero. It may be viewed as a function /o(tt) on 5' 2d_1 
that, besides, is S'0(d)-invariant. We may then rewrite the definition of A as 

T 



A = 

T— > 



lim^^jn fo(7r)Pt(Tro;dnj)dt. (9.2) 



Now, the existence of the limit follows from the fact that the Cesaro means of the transition probabilities 
Pt(7To; diz) tend in weak topology to the unique invariant probability measure ^(diz). Hence 

A = j /o(7r)/i(d7r) (9.3) 

and is independent of pq. The crucial input that allows to make the latter formula more explicit is the 
formula 

Lln^ = = £. (9.4) 



T T 



It implies that 



/o(tt) = Lln^T^2 = _ + Lh i, 

t 2 V p z 

x 1 fln(l + x 2 ) for d = l 

t + 2 | In (1 + x 2 + y 2 ) for d > 2 ( ' ' 

in terms of the SO(d) invariants with L given by explicit formulae (|7.8|) . 1 5[) or (|7.24[) . If the functions 
ln(l + a; 2 ) in <i = 1 and In (1 + x 2 + y 2 ) in d = 2 that are homogeneous of degree zero on M. 2d \ {0} were 
smooth, then their contributions to the expectation with respect to the invariant measure on the right hand 
side of (|9.2|) would drop out by the integration by parts. The problem is, however, the lack of smoothness of 
those functions at p — and a more subtle argument is required. 



9.1. d = 1 case. In one dimension, Eq. (|9.3j) reduces to the identity 



A 



OO 

/ (r + \ LXXY ^ 1 + x2 )) r l( x ) dx ( 9 - 6 ) 



with r)(x) given by Eq. (|7.6[) . Since the latter integral represents the integration of a smooth function against 
a smooth measure on S , it converges absolutely. Consequently, the formula for A may be rewritten in the 
form: 

n 

A= lim [ (- + lLln(l + x 2 ))r](x)dx. (9.7) 

n— foo J \T 2 / 

— n 

Now the integration by parts and the formula L'r) = for the formal adjoint ]J defined with respect to the 
Lebesgue measure dx show that the term with L In Vl + x 1 drops out (for the cancellation of the boundary 
terms it is crucial that the integral is over a symmetric finite interval [— n,n]). We obtain this way the 
identity 

n oo 
1 • f 1 f 

A = — lim / xi](x)dx = — p.v. / xrj(x)dx (9-8) 

T n->oo J T J 

— n — oo 

where "p.v." stands for "principal value". The result may be expressed [21] by the Airy functions [T3"] : 

A = -^- + -^4- ln(Ai 2 (c)+Bi 2 (c)) for c= -. (9.9) 

2r 4rc^ dc v w; (4r(A + 2B))f V ; 

The number A + ^r is the Lvapunov exponent for the one-dimensional Anderson problem (|3.5p . recall relation 
(|3.1[) . It is always positive reflecting the permanent localization in one dimension. On the other hand, A itself 
changes sign as a function of r and A + 2i3 signaling a phase transition in the one-dimensional advection 
of inertial particles [33 . 
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9.2. d > 2 case. In two or more dimensions, Eq. (|9.3|) becomes 

A = j + ^Lln(l + x 2 +y 2 ])r 1 (x,y)dxdy, (9.10) 

where n(x,y) is the density of the invariant measure from Eqs. (|7.9|) or (|7.17jl . The asymptotic behavior of 
n{x, y) was established in Sec. 17.21 and Sec. 17.31 We show in Appendix [B] that it guarantees that the term 
with L In (1 + x 2 + y 2 ) may, indeed, be dropped from the expectation on the right hand side of Eq. (19.21) so 
that 

A = — J x n(x,y) dx dy , (9-11) 

where the integral converges absolutely as follows from the estimates 14[) and (|7.23l) . In general, there 
is no closed analytic expressions for the right hand side, unlike in the one-dimensional case. The results of 
numerical simulations for A, indicating its qualitative dependence on the parameters of the model, together 
with analytic arguments about its behavior when At — > oo with A/B = const, or when r — > with 
A/t = const, and B/t = const, may be found in [2ll4| rT7lf23| f24] . 



10. CONCLUSIONS 

We have studied rigorously a simple stochastic differential equation (SDE) used to model the pair dispersion 
of close inertial particles moving in a moderately turbulent flow 2,17,23,24,29,33]. We have established the 
smoothness of the transition probabilities and the irreducibility of the dispersion process using Hormander 
criteria for hypoellipticity and control theory. For the projectivized version of the dispersion process, these 
results implied the existence of the unique invariant probability measure with smooth positive density as well 
as exponential mixing. The latter properties permitted to substantiate the formulae for the top Lyapunov 
exponents for the inertial particles used in the physical literature. In two space dimensions, we also showed 
that the complex-projectivized version of the dispersion process is non-explosive when described in the in- 
homogeneous variable of the complex projective space, unlike the real-pro jective version of the dispersion in 
one space dimension. A similar result was established in d > 3 for the complex-valued process built from the 
SO(d) invariants of the projectivized dispersion that was shown to stay for all times in the upper half-plane. 
These non-explosive behaviors are the reason why the numerical simulation of the processes in the complex 
(half-)plane could lead to reliable numerical results [2]. There are other questions about the models studied 
here that may be amenable to rigorous analysis. Let us list some of them: What about the expressions for 
the other 2d — 1 Lyapunov exponents (the 2d exponents have to sum to —d/r [TU])? Can one establish 
the existence of the large deviation regime for the finite-time Lyapunov exponents (the corresponding rate 
function for the top exponent was numerically studied in d = 2 model in [US]; it gives access to more 
subtle information about the clustering of inertial particles than the top Lyapunov exponent itself)? Is the 
SDE modeling the inertial particle dispersion in fully developed turbulence, that was introduced and studied 
numerically in OH], amenable to rigorous analysis? All those open problems are left for a future study. 



Appendix A. 

We establish here the expressions (|7.18[) and (|7.19j) for the SO(2d) -invariant normalized volume measure 
on 5 2d_1 for d > 3 (the same proof works also in d = 2, although there, the corresponding formulae are 
straightforward and well known). Note that for S 2 ^ 1 identified with the set of (p,x) G ^ 2d such that 
P 2 + X 2 = R 2 f° r fixed R, we may write 

dpo = const. 6(R — \/ p 2 + X 2 ) dpdx > (A-l) 

provided that we identify functions on S 2 ^ 1 with homogeneous function of degree zero on M. 2d \ {0}. Let 
us parametrize: 

p = O(p,0,...,0), X = 0(px,py, 0,...,0) (A.2) 

where p = \p\, x = \x\> x an d y are the .SO (^-invariants of Eq. (17.41) . and O € SO(d). Note that O^ 1 
is the rotation that aligns p with first positive half-axis of R d and brings x m t° the half-plane spanned by 
the first axis and the second positive half-axis, as required in Sec. 17.31 O and OO' , for O' rotating in the 
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subspace orthogonal to the first two axes, give the same (p, %). Let A be a d x d antisymmetric matrix, 
Aij = —Aji. Setting O = e A and differentiating Eqs. (|A.2j) at A = 0, we obtain: 

dp = (dp, pdA 21 , , . , , pdA dl ) , (A. 3) 

d X = (-pydA 21 ,pxdA 21 ,pxdA 31 + pydA 32 ,...,pxdA dl + pydA d2 ) 

+ (xdp + pdx,ydp + pdy, 0, . . . ,0) . (A. 4) 

Hence for the volume element, 

dpdx = p 2d ~ 1 y d ~ 2 dpdxdydA 21 ---dA dl dA 32 ---dA d2 . (A.5) 

The product of dAij gives, modulo normalization, the 50((i)-invariant volume element d[0] of the homo- 
geneous space SO(d) / SO(d — 2) at point [1]. Using the .SO^-invariance, we infer that 

dpdx = const. /? 2d_ V _5 'dpdx dy d[0] . (A.6) 
Substituting the last expression to Eq. (|A.1[) and performing the integral 

J 5(R - vV + x 2 ) P 2d ~ X d P = [S{R- P ^l + x 2 +y 2 ) P 2d ^d P 

r>2d-l 

(A.7) 



(1 + x 2 + y 2 ) d 

that collects the entire p-dependence in the integration against homogeneous function of zero degree, we 
obtain Eq. (|T. 19[) . modulo a constant factor that is fixed by normalizing of the resulting measure. 



Appendix B. 

We show here that 

J { L 9){x,y)r](x,y)dxdy = (B.i) 

in two or more dimensions, where g(x,y) = ln(l + .T 2 +j/ 2 ) and r](x,y) is the density of the invariant measure 
as defined by Eqs. (|7.9[) and (|7.17p . As mentioned in Sec.[5j the identity (|B.1[) does not follow immediately 
by integration by parts since the function g(x,y), is not smooth on S 2 ^ 1 . We shall then replace g(x,y) by 
the functions 

x 2 + y 2 \ _ / p 2 + x 2 
l + e(x 2 +y 2 )) ~ n Vp 2 + ex 2 

that are smooth on 5 2d_1 for e > 0. The identity (|B.lj) will follow if we show that 



g £ (x,y) = In ( 1 + -JL+VJ) = i n ( p +x ) ( B .2) 



Note that 



[Lg){x,y)ri{x,y)dxdy = lim (Lg e )(x,y) r)(x,y) dx dy . (B.3) 

e V J 



d y ge(x,y) 



(1 + (1 + e)(x 2 + y 2 ))(l + e{x 2 + y 2 )) ' 
2y 

(f + {1 + e)(x 2 + y 2 ))(l + e(x 2 + y 2 )) ' 
1 + (i + 2e)(y 2 - x 2 ) + e(l + e)(y 2 - 3x 2 ) 



so that 



d * 9e ( X ' V ) 2 (l + (l +e )( a; 2 +y 2 )) 2 (1+e(a .2 +y 2 )) 2 • 

d 2 o(xy) = 2 l + (l + 2 e )(^-y 2 ) + £ (l + e )(, 2 -3 y 2 ) 

o y g c [x,y) (1 + (1 + e)(a ,2 +y2)) 2 (1 + e(a .2 + y 2 )) 2 

1^(^)1 < - , 5 , |^<?e(^)l < 2W 



l + x 2 +y 2 ' 1 yaev ' yyl " l + a; 2 + y 2 ' 

\d 2 x g e (x,y)\ < fO, \d*g £ (x,y)\ < 10. (B.5) 

Using the explicit forms 1 5[) and (|7.24l) of the generator L, we infer that 

\(Lg e )(x,y)\ < C(l + \x\) (B.6) 
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with an e-independent constant C. Since the integral 

J(l + \x\) V (x,y)dxdy (B.7) 

converges due to the estimates (|7.14[) and (|7.23l) . and point- wise 

]im (Lg e )(x,y) = (Lg)(x,y) , (B.8) 

relation (IB. 31) follows from the Dominated Convergence Theorem. 
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